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We consider a random matrix model which describes the competition between chiral symmetry 
breaking and the formation of quark Cooper pairs in QCD at finite density. We study the evolution 
of the phase structure in temperature and chemical potential with variations of the strength of 
the interaction in the quark-quark channel and demonstrate that the phase diagram can realize a 
total of six different topologies. A vector interaction representing single-gluon exchange reproduces 
a topology commonly encountered in previous QCD models, in which a low-density chiral broken 
phase is separated from a high-density diquark phase by a first-order line. The other five topologies 
either do not possess a diquark phase or display a new phase and new critical points. Since these five 
cases require large variations of the coupling constants away from the values expected for a vector 
interaction, we conclude that the phase diagram of finite density QCD has the topology suggested 
by single-gluon exchange and that this topology is robust. 
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A number of early and recent model studies of finite density QCD jj]-|4j have suggested that quark Cooper pairs may 
form above some critical density and lead to 'color superconducting' matter. Although perturbation theory performed 
on single-gluon exchange suggests pairing gaps of a few MeV @, some recent calculations including non-perturbative 
interactions, either in the form of the Nambu Jona-Lasinio (NJL) model HH || or those induced by instantons j|,U, 
indicate gaps as large as 100 MeV. Theses values imply that color superconductivity may be relevant to the physics 
of heavy-ion collisions and neutron stars. 

Quark pairing singles out one color direction and thus spontaneously breaks ST/(3) co i or to SU(2) co i OI . The pattern 
of symmetry breaking may however be richer, since the formation of condensates in one (qq) channel competes with 
the breaking of chiral symmetry in the orthogonal (qq) channel. In an earlier paper || we formulated random matrix 
models for both chiral and diquark condensations in the limit of two quarks flavors and zero chemical potential. Our 
aim was to understand the phase structures which result from the competition between the two forms of order solely 
on the basis of the underlying symmetries. In this spirit, we constructed random matrix interactions for which the 
single quark Hamiltonian satisfies two basic requirements: (1) its block structure reflects the color SU(N C ) and chiral 
SU(2)l x SU(2)-r symmetries of QCD and (2) the single quark Hamiltonian is Hermitian. This last condition ensures 
the existence of well-defined relationships between the order parameters and the spectral properties of the interactions. 
In particular, condition (2) is obeyed by single-gluon exchange, which is generally regarded as the relevant description 
of QCD. Aside from conditions (1) and (2), the dynamics of the interactions does not contain any particular structure. 
In order to solve the model exactly, we described this "dynamics" by independent Gaussian distributions of matrix 
elements. 

In practice, what distinguishes the interactions from one another is their respective coupling constants, A and B, 
in the (qq) and (qq) channels. The ratio B/A measures the balance between chiral and diquark condensation forces. 
The absolute magnitudes of A and B play a secondary role. They introduce a scale in the condensation fields but do 
not affect the phase structure. We have shown in || that the condition (2) of Hermiticity forces B/A to be smaller or 
equal to N c /2. This constraint results in the absence of a stable diquark phase in the limit of zero chemical potential. 

The purpose of this paper is to extend our previous analysis to non-zero chemical potentials, for which the inter- 
actions cease to be Hermitian. Non-Hermitian interactions lead to considerable difficulties in numerical calculations 
in both lattice and random matrix theories. Standard Monte Carlo techniques fail because the fermion determinant 
in the action is complex and the sampling weights are no longer positive definite. Obtaining reliable results thus 
requires a proper treatment of cancellations between a large number of terms; a problem which has not yet found a 
satisfactory solution [p~p| — P~5f] - The present random matrix models possess exact solutions. The saddle-point methods 
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used in M to derive the free-energy as a function of the condensation fields remain valid at finite densities and give 
exact results in the thermodynamic limit of matrices of infinite dimensions. Here, we will use these methods to 
calculate the thermodynamics quantities as a function of the condensation fields and deduce the phase diagrams from 
the field configurations which maximize the pressure. 

We will discover that the pressure function has a simple analytic dependence which leads to polynomial gap 
equations. This situation reminds us of the gap equations obtained in a Landau-Ginzburg theory near criticality, and 
the random matrix approach is analogous in many respects. Both theories associate a change of symmetry with a 
change of state in the system. They are both mean-field and describe the dynamics of a reduced number of degrees 
of freedom N, where N scales with the volume of the system. In random matrix models, these degrees of freedom 
can be related to the low-lying quark excitation modes in the gluon background, i.e. the zero modes in an instanton 
approximation of that background Q In the Landau-Ginzburg formulation, the degrees of freedom correspond to 
the long-wavelength modes which remain after coarse-graining. There is, however, an essential difference in the 
construction of the two theories. A Landau-Ginzburg theory starts with the specification of an effective potential for 
the relevant degrees of freedom. A random matrix model starts at a more microscopic level with the construction of 
an interaction which, once integrated over, produces an effective potential. The integration can carry along dynamical 
constraints and thus restricts the allowed range of coupling constants. An example of such restriction is the Hermiticity 
condition in ||, which implies that B/A < N c /2. This condition is characteristic of the dynamics of the interactions 
which we consider and remains true at finite chemical potential. We will therefore take it into account in the following. 

In this paper, we consider random matrix models in which both chiral and color condensations can take place and 
study the resulting phase diagrams in temperature and chemical potential. We review the model of Ref. Q and 
discuss the general form of the pressure as a function of the condensation fields in Sec. H. We then solve the gap 
equations and analyze the six topologies that the phase diagram can assume in Sec. HI. We compare our results with 



those of QCD effective models, discuss the possibility of other symmetry breaking patterns and the extension to the 
case of non-zero quark current masses in Sec. [V. Section ^ presents our conclusions. 



II. FORMULATION OF THE MODEL 



A. The partition function 



The generalization of the matrix models introduced in Ref. to finite quark densities is straightforward. We 
represent the quark fields for each of the two flavors by ip 1 and ip 2 . The partition function at temperature T and 
chemical potential fi is then 



ip\ \ T ( H + (ttT + ifi) 70 + im 



x exp 



tiPa 

-U T + (ttT - ifi) 7 J - i 




(1) 



Here, TL is a matrix of dimension 4 x N c x N which represents the interaction of a single quark with a gluon background. 
Its measure, T>H, will be discussed below. The parameters m and 77 select a particular direction for chiral and color 
symmetry breaking and are to be taken to zero in the appropriate order at the end of the calculations. The current 
quark mass m is to be associated with the chiral order parameters (V'jV'i) an d ^2} ■ The complex parameter rj 
is to be associated with the order parameter for Cooper pairing, (i/>J '-Pa^i) 1 m which Pa = iC^\2 (C is the charge 
conjugation operator) selects the quark-quark combinations which are antisymmetric in spin and color, i.e. a chiral 
isosinglet, Lorentz scalar, and color 3 state. Note that in order to permit the construction of correlations between the 
fields ipi and 1P21 we have transposed the single quark propagator in the second flavor, hence the upperscript T. 

The interaction Ti. is intended to mimic the effects of gluon fields and thus explicitly includes the desired chiral and 
color symmetries. Of central interest is single-gluon exchange, which has the chiral block-structure 



'This relation is explicit in the random matrix models where only chiral symmetry is considered |l6|-[l8t; it is expected to 
remain true once color symmetry is also taken into account M. 
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H -( ° W ™\ (2) 
where W sge has the spin and color block-structure of a vector interaction, 

4 8 

w sge = j2J2 a i® xa ® A,M *- (3) 

/i=l a=l 

Here, cr+ = (l,^)^ are the 2x2 spin matrices, and A a denote the N c x N c Gell-Mann matrices. The A^ a are real 
N x N matrices which represent the gluon fields. 

Since we want to explore the evolution of the phase structure as the balance between chiral and diquark conden- 
sations is changed, we consider the larger class of Hcrmitian interactions to which single-gluon exchange belongs. As 
noted in the introduction, this choice is motivated by the fact that Hermitian interactions have a clear relationship 
between the order parameters and the spectral properties, in the form of Banks-Casher formulae We write an 
Hcrmitian interaction H. as an expansion into a direct product of the sixteen Dirac matrices Tc times the N£ color 
matrices. The matrix elements are given by 

16 N 2 C 

HxkxkiKjfil = ( r c) Xi; K j Yj A «/5 ( A \Z) kl > (4) 
C=l a=l 

where the indices (A,k), and (a, (3) respectively denote chiral, spin, and color quantum numbers, while (k,l) 

are matrix indices running from 1 to N. The A a represent the color matrices A a when a < — 1 and the diagonal 
matrix (S c ) a p = S a p when a = N%. The normalization for color matrices is Tr[A a A & ] = 28 a b and Tr[5^] = N c ; the 
normalization of the Dirac matrices is Tr[TcTc*/] = 45cc- 
The random matrices 

A Ca 

are real when C is vector or axial vector (C = V, A) and real symmetric when C is 
scalar, pseudoscalar, or tensor (C — S,P,T) B. Their measure is 



Ca Xk 



Ca Xk 



(5) 



where T>A^ are Haar measures. Here, (3c = 1 for C — V,A and (3c = 1/2 for C = S,P,T. We want to mimic 
interactions which in a four-dimensional field theory would respect color SU (N c ) and Lorentz invariance in the vacuum. 
Therefore, we choose a single variance Sc a for all channels which transform equally under color and space rotations. 

The temperature and chemical potential enter the model in Eq. (^) through the inclusion of the first Matsubara 
frequency in the single quark propagator. Such T and fi dependence is certainly oversimplified but none the less 
sufficient to produce the desired physics. Our purpose is to understand the general topology of the phase diagram 
and not to provide explicit numbers. More refined treatments including, for instance, all Matsubara frequencies 
would modify the details of the phase diagram and map every (/i, T) coordinate to a new one. However, any such 
mapping will necessarily be monotonic and will conserve the topology. We note in particular that we do not assume a 
temperature dependence of the variances, i.e. we neglect the T-dependence of the gluon background. A non-analytic 
behavior in T should only arise in the contribution to the thermodynamics from the degrees of freedom related to 
chiral and diquark condensations. We thus expect a realistic T-dependence of the gluon background (and of the 
variances) to be smooth and not to affect the overall phase topology. 

We choose the signs of the T- and /i-dependencies to mimic a diquark condensate which is uniform in time, i.e. 
which does not contain a proper pairing frequency. In a microscopic theory formulated in four-momentum space, the 
absence of a proper frequency leads to pairing between particle and hole excitations with energies which are symmetric 
around the Fermi surface. We simulate this effect by selecting opposite T-dependences for the fields tp 1 and tjj^ , while 
maintaining the same /^-dependence. 

B. The pressure function 

The integration over the random matrix interactions is Gaussian and can thus be performed exactly. Following the 
procedure of Ref. ||, we use a Hubbard-Stratonovitch transformation to introduce two auxiliary variables a and A, 
to be associated with the chiral and pairing order parameters respectively 
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An integration over the fermion fields then reduces the partition function to 



Z(fi,T) = J dcrdAexp [— 4Ml(cr, A)] 



(6) 



(7) 



where f2(<r, A) is the negative of the pressure, P(<r, A), per degree of freedom and per unit spin and flavor, 
f2(<7, A) = -P(cr, A) 

= AA 2 + Bo 2 - i|(JV c - 2) log [((a + to - M ) 2 + T 2 ) ((a + to + i u) 2 + T 2 ) 

+2 log [((o- + to - ^) 2 + T 2 + | A + r?| 2 ) ((a + m + + T 2 + |A + r?| 2 )] J, 



(8) 



where we have dropped the prefactors 7r in the temperature dependence for simplicity. Here, the coupling constants A 
and B are weighted averages of the Fierz coefficients /^ and f£ a obtained respectively by projecting the interaction 
Fc ® A a onto chiral and diquark channels, 



\Ca ) \Ca J 



(9) 



To make contact to microscopic theories, we note that a small coupling limit in either channel corresponds to a small 
Fierz constant and hence to large parameters A or B. This limit favors small fields A or a. Because we can always 
rescale the condensation fields by either y/~A or y/B, the only independent parameter in Eq. (^) is the ratio of B/A, 
which by virtue of Eq. (j^) is a measure of the balance between the condensation forces. Again, Hermitian matrices 
H satisfy B/A < N c /2. 

The mass m and the parameter r\ explicitly break chiral and color symmetries. They act as external fields which 
select a particular direction for the condensation pattern, and should be taken to zero at the end of the calculations. 
(They can also be kept constant to study the effect of a small external field, a point which we take in Sec. IV.) They 



are useful for obtaining the order parameters from derivatives of the partition function. In the thermodynamic limit 
N -> oo, Z(n, T) in Eq. (g) obeys 



lim log Z [fx, T) 

N^oo 



N 



lim 4N min CT! A {0(cr, A)} 



(10) 



where the right side represents the global minimum of f2 in Eq. 
and diquark condensations are given as 



for fixed fi and T. The order parameters for chiral 



= lim lim 



1 dlogZ 



m->ojv^oc ANfN dm 



(il> T '4') = ~ li m 1™ 



f d log Z 



?7,?7*— 



Tj,r;*^0 JV^oo 47V drj* 



(11) 



m=0 



where the number of flavor is Nr — 2. Note that the thermodynamic limit N — > oo must be taken first before the 
small field limit to, 77 — ► 0, see pl|]. Given the to- and 77-dependences of the log terms in Eq. (ph, we have 



(12) 



where cr(n, T) and A(/i, T) are the condensation fields which minimize Q(cr, A) for fixed \x and T. 

Equation (pj) is the main expression from which we will deduce the phase structure as a function of the ratio B/A. 
Its form is very simple to understand. The quadratic terms correspond to the energy cost for creating static field 
configurations with finite a and A. The log terms represent the energy of interaction between the condensation fields 
and the quark degrees of freedom. They can be written in a compact form as Tr[log S(a, A)] = logdet S(a, A), where 
Tr is a trace in flavor, spin, and color, and S is the single quark propagator in a background of a and A fields. 
Substituting the Mastubara frequency T — > ipi, it becomes clear that the poles of S in correspond to the excitation 
energies of the system. From Eq. (j§), we see that two colors develop gapped excitations with 



Pi 



= ±V(^tm) 2 + a 2 , 



(13) 
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where the plus and minus signs respectively correspond to particle and antiparticle modes. The N c — 2 remaining 
colors have ungapped excitation with 

P4 = ±|ffT/i|. (14) 

It is important to recognize that the potential Q in Eq. (^) contains the contribution to the thermodynamics from 
only the low energy modes of Eqs. and (fl4|). The right side of Eq. (|J) thus corresponds to the non-analytic piece 
in the thermodynamic potential which describes the critical physics related to chiral and color symmetry breaking. 
In a microscopic model, the right side of Eq. (8) would also contain a smooth analytic component Q reg (p,,T) which 
arises from all other, non-critical degrees of freedom in the system. Our model does not contain this contribution and 
thus cannot be taken as a quantitative description of bulk thermodynamics properties; the model is constructed to 
describe the critical properties and should be used as such. 

We show in the next section that the forms of the potential f2(cr, A) in Eq. (^) and of the associated excitation 
energies in Eqs. ([l3| ) and ( |l4|) are sufficient to produce a rich variety of phase diagrams which illustrate in a clear way 
the interplay between chiral and color symmetries. In particular, we will find that single-gluon exchange reproduces 
the topology obtained in many microscopic models of finite density QCD, see for instance @-@|. 

III. EXPLORING THE PHASE DIAGRAMS 

We noted earlier that the potential Q in Eq. (||) is very similar to a Landau- Ginzburg functional. Although il(er, A) 
is not an algebraic function, it is equivalent to a polynomial of order a 6 and A 4 : The gap equations reduce to coupled 
polynomial equations of fifth order in a and third order in A. Generically, four types of solutions exist: 

• (i) the 0-phase, the trivial phase in which both a and A vanish, 

• (ii) the x-phase, in which chiral symmetry is spontaneously broken but A = 0, 

• (iii) the A-phase, in which color symmetry is spontaneously broken but a = 0, 

• (iv) the %A phase, a mixed broken symmetry phase in which both fields are non- vanishing, a ^ and A ^ 0. 

The ^A-phase is thermodynamically distinct from the \~ an d A-phases and is not a "mixture" of these two phases. 

At a given /i, T, and B/A, each of the solutions (i) - (iv) can either be a minimum or a saddle-point of f2. The 
complex flow of these solutions with the variation of B/A leads to large variety of phase diagrams. The phase 
structures can however be grouped according to their topologies as shown in Figs, [j]- [| Figure || shows the case 
of smallest values B/A, which favor chiral over diquark condensation. B/A then increases continuously from Fig. |^ 
to ^. We now discuss Figs. 1-6 in the six following subsections. We indicate in parentheses the corresponding ranges 
of B/A and discuss their limits in the text. 

Pure chiral condensation: (0 < B/A < 0.1397V C r < B/A < 0.418 for N c = 3; see Fig. 1). Generically, 
the global minimum of f2 in Eq. (|^) is realized by field configurations which maximize the log terms while keeping 
reasonably low values of the quadratic terms A A 2 + Bo 1 . We consider first the limit A ^> B. The minimum of 
must then correspond to A = in order to avoid the large energy penalty ~ AA 2 . No diquark condensation occurs 
in this case, and color enters only as a prefactor N c in the number of degrees of freedom. This becomes clear if we 
absorb N c into B and set B to 1; we then recover the potential SI studied in the chiral random matrix models pcj ] 
which neglect color altogether. 

We briefly recall the phase structure in this case. The gap equation for the chiral field, d£l(a, Q)/da = 0, has a 
trivial root a = and four other roots which satisfy the following quadratic equation for a 2 

N c {pi 2 - T 2 ) + B(fi 2 + T 2 ) 2 + (2B(T 2 - - N c )a 2 + Bcr 4 = 0. (15) 

In the high temperature limit T 3> /i, both roots in a 2 are negative. The only real solution of the gap equation is 
(7 = 0, and the system is in the symmetric phase. Decreasing T for fixed /i, one encounters a line of second-order 
phase transitions 

L Xt 2 = N c ([j, 2 — T 2 ) + B(fj? + T 2 ) 2 = 0, (16) 



2 We have rescaled in each figure the units of T and /i for clarity. 
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where one of the roots a 2 in Eq. ([15]) vanishes. Below £ x ,2, the trivial root becomes a local maximum of 0(<r, 0), and 
we have a pair of local minima at the real roots a = ±<toj where 




^o = Ui-J"+^+ v c „; • (17) 



These roots correspond to a chiral broken phase. They become degenerate with a = on the second-order line i x ,2, 
where the potential O scales as 0(cr, 0) — O(0, 0) « it 4 . Thus, the critical exponents near L Xj 2 are those of a mean-field 
</> 4 theory. 

The second-order line ends at a tricritical point (^3, T3) at which all five roots of the gap equation vanish and where 



0(<r, 0) — O(0,0) ss cr 6 , giving now critical exponents of a mean-field </> 6 theory. From Eq. Jl5|), this happens when 
2B(Tf - (if) - N c = 0, which with the use of Eq. © gives 





K = \hhVv2-\, t 3 = \/ 7 ^V^ 2 + 1 - ( 18 ) 



For fi > fis, the transition between the chiral and trivial phases is first-order and takes place along the line of equal 
pressure 




£x,i^|l + Jl-16M 2 T 2 ^-log 




N 2 

+N C log [fi 2 + T 2 ] + B(n 2 - T 2 ) = 0. (19) 
This line intercepts the T — axis at /i = p,\ which obeys 

1 + f + «f> = ' < 2 »> 

This gives m = 0.528 s/N c /B, or = 0.914/-/B for iV c = 3. 

L x .i in Eq. ( |H| ) is a triple line. To see this and clarify the character of (/13, T3), it is useful to consider the effect of 
a non-zero quark current mass to p0[ . A mass m/0 selects a particular direction for chiral condensation. If we now 
consider the three-dimensional parameter space (n, T, to), the region delimited by L x ^ and L Xi i in the plane to = 
thus appears to be a surface of coexistence of the two ordered phases with the chiral fields ±00 of Eq. (|l^) . Along 
L x ,i, this surface meets two other 'wing' surfaces which extend symmetrically into the regions to > and to < 0. 
Each of the wings is a coexistence surface between one of the ordered phases whose chiral field continues to a = ±00 
as to — ► and the high temperature phase. Hence, L x 1 marks the coexistence of three phases and is a triple line. 
The three phases become identical at (^3,13), which is thus a tricritical point. The second- and first-order lines L x .i 
and L x ,2 join tangentially at (fi 3 ,T 3 ), see p2[ . 

The onset of diquark condensation modifies the topology that we have just described. This takes place for coupling 
ratios B/A > 0.139A C to which we now turn. 

The QCD case: (0.139A C < B/A < a x {N c ) or 0.418 < B/A < 1.05 for N c = 3; see Fig. 2). To understand the 
conditions for the onset of diquark condensation, consider a pure diquark phase by setting cr = in Eq. (j^) . The gap 
equation dft/dA = then has three solutions: A = 0, and A = ±Ao where Ao = \/2/A — — T 2 . In the high T 
and n phase, only the trivial root is real and the system is in the symmetric phase. Inside the semicircle 

L A:2 = 2/ A - - T 2 = 0, (21) 

the trivial root becomes a maximum of O(0, A), and we have a pair of two real minima A = ±A for which color 
symmetry is spontaneously broken. The roots ±Ao go continuously to zero as one approaches £a,2, which is thus a 
second-order line provided no other phase develops. The fate of other ordering forms depends on the pressure in the 
diquark phase, 

P A = -0 A = A(/i 2 + T 2 ) - 2 + (JV C - 2) log (/i 2 + T 2 ) + 2 log f | \ . (22) 

The maximum pressure Pa is reached on La, 2, where it also equals the pressure of the trivial phase. The condition 
for the onset of diquark condensation is then clear: The semicircle must lie in part outside the region occupied by 
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the chiral phase. Then, the maximum pressure in the diquark phase is higher than that of the chiral phase, and the 
diquark phase is stable near the semicircle. 

To find the minimum ratio B/A for which this condition holds, we compare the dimensions of the chiral phase to 
the radius of the semicircle, /i scm i = \/2/A. When B/A -C 1, the line La, 2 lies well inside the boundaries of the 
chiral phase, whose linear dimensions are T = <J N c /B (along /i = 0) and fix = 0.528 ^JNJB (along T = 0). The 
onset of diquark condensation requires therefore B/A to be large enough that the semicircle crosses the first-order 
line between chiral and trivial phases. This takes place along T = when /x sem j = Hi, or 



y/2/A = 0.528^/^/5, (23) 

which gives B/A = 0.1397V C , or B/A = 0.418 for N c = 3. 

For larger values of B/A, the diquark phase exists in a region delimited by La, 2, on which it coexists with the 
symmetric phase, and by a first-order line of equal pressure with the chiral phase. This phase diagram is realized for 
ratios including that corresponding to single-gluon exchange, B/A = N C /(2(N C — 1)) (or 3/4 for N c = 3), which is the 
ratio taken in Fig. |^. We have verified that the two segments of first-order lines, between the chiral and the diquark 
phases on the one hand and the chiral and trivial phases on the other hand, join tangentially. Figure |?j shows the 
chiral and diquark fields as a function of [i and T. It is worth noting that the chiral field vanishes with a square root 
law near the second-order line L Xi 2 and is discontinuous along the first-order line L Xi %. The diquark field vanishes 
with a square root law all along the semicircle line. 

The topology in Fig. ^ can also be summarized by a simple counting argument. Consider the difference between 
the pressures in the diquark and the chiral phases along the T = axis, 

AO = n x - fl A = N c - 2 + (A + B)[± 2 - 2 log + N c log (24) 

This difference varies from Afl = N c (l + (2B)/(N C A) + log[(2B)/(iV c A)]) as fi -> y/2/A, to AQ = (N c - 2)log^ 2 
as fi — > 0. We have just argued that Afi can be negative in the former limit if B/A > 0.139^. In the latter limit, 
however, Af2 is always negative for N c > 2 and the system is necessarily in the chiral phase. That the diquark 
phase must be metastable for small /1 is clearly a consequence of the fact that chiral condensation uses all N c colors 
while diquark condensation uses only colors 1 and 2. This counting argument was mentioned in early models of color 
superconductivity || and is here completely manifest. For N c > 2, the diquark phase must appear at densities higher 
than those appropriate for the chiral phase. 

Onset of the mixed broken symmetry phase: (ai(N c ) < B/A < N c /^/S or 1.05 < B/A < 1.06 for N c = 3; 
see Fig. 3). The chiral solution ctq, Eq. (p7|), ceases to be a local minimum of 0(cr, 0) for sufficiently large /z. It turns 
into a saddle point along a line where the second derivative d 2 Vl/dA 2 vanishes, 



{N C A - 2B) N c + y/N 2 - 1QB 2 ii 2 T 2 - 8B 2 ^ 2 = 0. (25) 



A new minimum of 0(<7, A) with both a ^ and A / develops in the region to the right of L m j X . This new 
minimum corresponds to a new phase, the xA-phase, which competes with the diquark phase. When the xA-phase 
first appears along L m j X , it has the same pressure as the chiral phase. Therefore, the condition for this new phase to 
realize the largest pressure is that the instability line L m i x lies in the region spanned by the chiral phase. Then, along 
L m j x , the new phase necessarily has a pressure that exceeds that in the diquark phase and is favored. 

To see what ratios B/A are needed for the mixed broken symmetry phase, consider the point where L m ; x meets the 
T = axis, 



2B V 2B 



When B/A <C 1, ^ m ix ~ \/N c /B/ y/4B/A 3> y/N c /B and L m i x lies well outside the chiral phase. L m i x then moves 
towards the x-phase as B/A increases. It first crosses the first-order line between chiral and diquark phases when 
Mmix = Mi; where [i\ is the point of equal pressure which obeys Ail = in Eq. (p4|). This condition gives a ratio 
B/A = a\(N c ), where ai is a non-trivial function of N c which we study in the appendix. Here, we note only that 
a 1 (3) = 1.05 for three colors. When B/A > ax(N c ), we find that the mixed broken symmetry phase develops in the 
wedge bordered by the second-order line L m ; x of Eq. (g5j) , and a first-order line on which it coexists with the diquark 
phase. The phase diagram is shown in Fig. ^, while an expanded view of the wedge of mixed broken symmetry and 
of the first-order line near (//3, T3) are respectively shown in Figs. || and ^. 
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A new critical point: (N C /VS < B/A < a 2 (N c ) or 1.06 < B/A < 1.163 for N c = 3; see Fig. 4). As B/A increases 
above ai(N c ), the semicircle 2/A = fj 2 + T 2 grows relative to the chiral phase, and eventually reaches the tricritical 
point (/i 3 ,T 3 ) ofEq. ( |l8| ) when B/A = N c /V8. At this stage, the segment of first-order line between chiral and trivial 
phases disappears. For B/A > N c /y/8, the semicircle meets the second-order line between chiral and trivial phases, 
L Xi 2 in Eq. (16]), at a new critical point 





T* = \H[1 + Tnh (27) 



which now separates three phases as shown in Fig. [|. 

Coexistence of four phases: (a 2 (N c ) < B/A < N c /2 or 1.163 < B/A < 1.5 for N c = 3; see Fig. 5). With still 
higher coupling ratios B/A, the wedge of mixed broken symmetry in Fig. ^ grows in size relative to the other phases. 
Its tip reaches the new critical point (114, T4) when B/A satisfies the condition 

B rAT ,_ 4N c -N? /2 V2N^4 

A =a2{Nc) = W^Q ' { ) 

the derivation of which we detail in the appendix. We just note here that «2(3) = 1.163. This coupling ratio also 
marks the appearance of two new critical points, as illustrated in Fig. ||. First, when B/A > a.2(N c ), the point 
(/X4, T4) characterizes the coexistence of all four phases, and it has thus become a tetracritical point. The pressure of 
the system at (^4, T4) has the form f2(cr, A) — O(0, 0) « aer 4 + 6A 4 + ca 2 A 2 , where a, b, and c are constants detailed 
in the appendix. We note in particular that none of the four second-order lines in Fig. 5 join tangentially at the 
tetracritical point. 

The second point appearing above B/A = a2(N c ) is a tricritical point (n3 m ,T^ m ) which lies on the boundary 
between the %A- and the A-phases. Its origin can be understood from the similarity between the characters of the 
mixed broken symmetry phase and the chiral phase. The diquark field in the xA-phase is 

( 1 + yi + 4AVa 2 2 2 2 V /2 cm\ 
A(<t) = I fj, —T — a J . (29) 



Substituting A in fi(cr, A) for this expression, we find that the gap equation for the chiral field, dfl(a, A(a))/da = 0, 
has five roots. Despite the fact that A(<r) ^ 0, the dynamics of these roots as a function of n and T is identical to 
that of the roots in the pure chiral condensation case discussed at the beginning of this section. In particular, the 
transition to the A-phase starts out second-order near (/i4,T4) and takes place along a line where three of the five 
roots vanish. This second-order line is thus determined by the condition 



d 2 Q , ., 



= 0, (30) 



which is the analog of Eq. ([TJ]) and gives 

L xA ^ A . 2 = (N c - 2)( M 2 - T 2 ) + (B- A(l - V)) (M 2 + T 2 ) 2 = 0. (31) 

In analogy with the pure chiral case, the transition becomes first order at high /1; second- and first-order segments 
join tangentially at a tricritical point (/i3m,T3m)- 

We determine the location of this point as follows. We observe that, as in the pure chiral broken case, the potential 
along the second-order line £ x a^a,2 nas the form fi(cr, A(ct)) — J7(0,0) w c 4 . (See Eq. (^0|).) This scaling form 
becomes f2(<7, A(cr)) — f2(0, 0) s» cr 6 at (/J,3 m , Tbrn) where all five roots vanish. Thus, the tricritical point is located on 
the line £ x a^a.2 at the point where 

0(*,A(<7)) =0, (32) 

which is the analog of the condition 2B(T 2 — n 2 ) — N c — for the pure chiral case. We have solved Eq. (^) numerically 
to determine (/Lt3 mj T3111) in Figs. ^ and |[ 

Disappearance of the chiral phase. (02 (-/V c ) < B/A; see Figs. 6 and 10). The wedge of mixed broken symmetry 
gains space in the (n,T) plane at the expense of the chiral phase as B/A increases above ct2(N c ). The chiral phase 
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eventually shrinks to a vertical line along the /i = axis when B/A — N c /2, as shown in Fig. [T(j. That the chiral 
phase is still present for \i — can be understood as follows. When B/A = N c /2 and /i = 0, both condensation fields 
appear in Eq. (^) in the combination a 2 + A 2 . A pure diquark solution (cr, A) = (0, Ao) can thus always be rotated 
into a pure chiral solution (cr, A) = (Ao, 0) ||, and thus does not represent an independent phase. This symmetry is 
however no longer present as soon as fj, ^ 0, in which case the diquark phase becomes thermodynamically independent 
from the chiral phase. 

Recall that the ratio B/A = N c /2 is the maximum ratio that our model can realize ||. It is instructive however to 
explore higher coupling ratios, although they do not necessarily describe physical situations, by setting B/A > N c /2 
by hand in Eq. (||). We find that higher B/A force the diquark phase to grow in size at the expense of the mixed 
broken symmetry phase. Figure || shows one example with B/A = 1.8 and N c — 3. The similarities between the 
critical properties of the mixed broken symmetry phase and those of the pure chiral broken phase are now clear: 
Compare for instance the critical lines between Figs, p] and @. 



IV. DISCUSSION 



Baryon density discontinuity. We have argued that the potential fl in Eq. (8) represents the non-analytic 
contribution to the thermodynamics which is directly associated with the breaking of chiral and color symmetry. The 
potential 17 should therefore not be used too literally in computations of the bulk properties of the phases we have 
encountered. It may however give reasonable estimates of discontinuities near a phase transition. For instance in the 
pure chiral case, the random matrix models of Ref. ||(]] estimate that the baryon density ns = — (l/3)<9f2(/Lt, T)/<9/z 
changes discontinuously at the first order point along T = by an amount Ans ~ 2.5 no. Here, uq = 0.17fm -3 is the 
density of normal nuclear matter. The result for Ans relies on an evaluation of the number of degrees of freedom N 
from instanton models and seems to be a reasonable estimate of the baryon density discontinuity [^of . 

Ans is modified by the presence of diquark condensation. We consider single-gluon exchange with N c ~ 3, which 
realizes B/A = 3/4, and work in the limit T = 0. Taking the derivative of Eq. ( pi| ) with respect to fi, we find the 
discontinuity in baryon density at the point [i\ of equal pressure between chiral and color phases to be 

A A - N x = + IbA . (33) 

Here, /ii obeys fl x = S1a in Eq. (|24|); we have /ii = 0.87 /\[B. Hence, — N x ~ 6Ay/B in the appropriate unit of 
inverse chemical potential. Were no diquark condensation to occur, we would have found a discontinuity 

where [i\ is now the point of equal pressure between chiral and trivial phases. We want to keep the same B as in 



Eq. (33) so as to compare two situations which realize the same vacuum chiral field Co = y/^/B. The condition of 
equal pressure in Eq. (|2(]) gives then [i\ = 0.914/y A B and we have Ao — N x ~ 8.4 */B. Thus, we find that diquark 
condensation reduces the discontinuity in baryon density by roughly twenty five percent. 

Away from the chiral limit. We now turn to study the effects of a small quark current mass m in Eq. (|sj) . For 
m 7^ 0, chiral symmetry is explicitely broken, and the chiral condensate (ip'tp) ceases to be a good order parameter. 
This affects the phase diagrams in Figs. 1-6 in a number of ways. We consider the effect of a small mass m chosen 
so that m ~ 10 MeV in units for which the vacuum chiral field is a w y/3/B ~ 400 MeV and illustrate a few cases in 
Figs. 11-14. 

Figure 11 shows the limit of small ratios B/A which favor chiral over diquark condensation. Since (ip^ip) 1S n0 longer 
a good order parameter, any two given points in the phase diagram can be connected by a trajectory along which no 
thermodynamic discontinuity occurs. It results that the second-order line £ x ,2 in Eq. ( pL6| ) is no longer present when 
m 7^ 0. There remains, however, a first-order line, which ends at a regular critical point (/x c ,T c ). This point can 
be located as follows. Along the first-order line, the potential f2(cr, 0) has two minima of equal depth separated by 
a single maximum. All three extrema become degenerate at the critical point (/x C) T c ), past which f2(<7, 0) possesses 
only one minimum. The location of (jj, c ,T c ) can thus be determined from the condition that f2(cr, 0) scales as 

n(tr, 0)- 17(0,0)- (a - a ) 4 , (35) 

at (/Li c , T c ) and for small deviations | cr — ctq | . The location of (/i c , T c ), as well as ao, are then determined by requiring 
the first three derivatives of Q (cr, 0) to vanish at the critical point and for cr = gq. A small mass tends to increase the 
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pressure in the low density 'phase' with respect to that in the high density 'phase'. Its results that a mass m displaces 
the first-order line £ Xj i of Eq. (|l9|) to higher [i by an amount linear in m. The pressure increase also delays the onset 
of diquark condensation; the ratio B/A for which the diquark appears first increases linearly with m. 

Figure 12 shows the case of QCD: there is now a second-order line La, 2 which separates the high temperature phase 
from a mixed broken symmetry phase; the dominant effect of m on the diquark phase in Fig. 2 is to produce a small 
chiral field a ~ m. The effect on the thermodynamics in both the chiral and the trivial phases is, however, second 
order in m, and so is the displacement of the second-order line La, 2 from Fig. 2 to Fig. 12. 

The evolution of the phase diagram for higher B/A parallels the evolution we have outlined for m = in Figs. 3-6. 
The phase with a finite diquark field grows until it eventually reaches the critical point (/i c , T c ). Higher ratios lead to 
a phase structure in which the first- and second-order lines intersect as shown in Fig. 13. A wedge of another mixed 
broken symmetry phase, initially with a large chiral field and a small diquark field, appears at still higher B/A on 
the left of the first-order line. This wedge grows in size until it encounters the boundary of the other mixed broken 
symmetry phase. The second-order lines then merge into a single continuous line as shown in Fig. 14. This line is the 
locus of point for which 



d 2 n{a, A) 



OA 2 



= 0, (36) 

A=0 



and on which a chiral solution with a vanishing diquark field turns into a saddle-point of f2(c, A). Inside this boundary, 
the global minimum of Q(cr, A) describes a single mixed broken symmetry phase. This phase again exhibits properties 
similar to those of a chiral phase. It contains in particular a first-order line which ends at a critical point (n C m-,T cm ) 
at which the potential fl scales as f2(<7, A(cr)) — £7(0, 0) ~ (a — <To) 4 , where A(cr) is the diquark field for fixed a. 

To summarize, the effects of a small mass is linear for the chiral and mixed broken symmetry phases and quadratic 
for the diquark and trivial phases. The chiral field no longer represents a good order parameter, and the second-order 
lines of vanishing second derivatives with respect to cr, i.e. L X 2 in Eq. ( [l6| ) and L m - lx in Eq. (EH), disappear. The 
tricritical points become regular critical points while the tetracritical point in Fig. 5 disappears. 

The N c = 2 and N c — > 00 limits. In order to make connection with known results and to obtain some insight 
on the dependence of the phase structure on the number of colors, it is interesting to consider the limits N c = 2 and 
N c — ► 00 . First of all when N c — 2, the Hermitian matrix models realize only a single coupling ratio B/A — 1 g|. We 
focus on single-gluon exchange for which the potential in Eq. (||) becomes 

n 2 (a, A) = A(a 2 + A 2 ) - Iog[(<r - [if + A 2 + T 2 } - log[(a + y) 2 + A 2 + T 2 }. (37) 

A rotational symmetry appears at /1 = as ^2 depends on chiral and diquark fields via the combination a 2 + A 2 . 
Diquark and chiral phases do not in this case represent independent states. We find for the combined fields that 

cr 2 + A 2 = j - T 2 , (38) 

below a critical temperature T c — y/2/A, while symmetry is restored above T c . This rotational symmetry is a 
consequence of the pseudo-reality of ST/(2)-QCD, a property by which the Dirac operator D = i Ylua J^^A^a + m 
commutes with r-iCy'K where T2 is the antisymmetric 2x2 color matrix and K the complex conjugation operator. 
For p, = 0, this property permits one to arrange color and flavor symmetries into a higher Si/ (4) symmetry, see for 
instance [p3|-p5|. 

At finite /z, however, the 5/7(4) symmetry is explicitly broken. The global minimum of f2(cr, A) always has a = 0, 
and the system prefers diquark condensation over chiral symmetry breaking. This results agrees with instanton 
models PJ26J , and lattice calculations J27J . We now have a second-order phase transition from an ordered state with 
A 2 = 2 /A — fj 2 — T 2 in the low T and /1 region to a symmetric phase at high T and fi. 

The opposite limit N c — * 00 is more subtle. In microscopic models, it is expected on general grounds that the 
quark-quark interaction is suppressed with respect to the qq channel by powers of 1 /N c |26| , |28|| . In the present model, 
we observe that diquark condensation disappears as N c —> 00 if the interaction is single-gluon exchange. Its coupling 
ratio B/A = N C /(2N C - 2) — > 1/2 as N c -> 00. Thus, we have B/A < 0.139iV c and the only possible topology 
for the phase diagram is that in Fig. 1. Hence, no diquark condensate forms. Other Hermitean interactions with 
B/A ~ 0(N C ) can however explore the full range of phase diagrams which display diquark condensation. The actual 
number of possible toplogies is however reduced to five, as a.2(N c ) — > N c /^/8 for N c — > 00 (see Appendix) and Fig. 4 
can no longer be realized. 

We summarize the variation with N c of the coupling ratios characterizing a change in topology and of the ratio 
realized by single-gluon exchange in Table 1. 
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Comparison with a microscopic model: In contrast to microscopic models, the random matrix interaction 
does not lead to a logarithmic instability of the gap equation near the Fermi surface. To clarify this effect, we 
consider diquark condensation in a chiral symmetric phase at T = and compare the random matrix approach to a 
microscopic model. We choose here the NJL study of Berges and Rajagopal M. In the random matrix formulation 
the gap equation 

2AA = (39) 
has a non-trivial root A ~ y/jjL c — \i with /Lt c = ^J2/ A. The model of |(| gives an equation of the form 

2A'A = A f°° dqq 2 F 4 (q) [ 1 + 1 ) , (40) 

where F(g) is an appropriate form factor which falls off over q ~ A ~ O(Aqcd). Clearly, the right side of Eq. ( |40| ) 
contains a singularity for q ~ fi and A — > which is absent in Eq. (^9|). This singularity has two consequences. First, 
the behavior of the diquark field for large fi depends sensitively on the form factor. As A — > 0, the singularity at q ~ fi 
in the right side of Eq. (EG) gives 2A' w fi 2 F 4 (fi) log A/x/A 2 to logarithmic order. Thus, instead of the square root 
behavior A(/i) ~ (Me - m) of the random matrix approach, A now vanishes as /i — * oo as A(/i) oc exp(— c / (/j 2 F 4 (fi))) , 
where c is a constant. f\ This tail is exponentially sensitive to form factors; with regards to establishing the general 
topology of the phase diagram, such behavior is not qualitatively different from A = 0. The second consequence of 
the logarithmic singularity in Eq. ( |40| ) is that A(/i) must be non-monotonic for intermediate (i. This profile can be 
seen by drawing in the plane (A, /i) the lines of constant height for the right side of Eq. (f40|): these lines must go 
around the point (/uo,0) where the logarithmic singularity is the strongest as A — > 0, and A(/i) reaches a maximum 
at [1q. 

Color dependency in the chiral fields. The patterns of symmetry breaking which we have described can 
become even richer if we allow the chiral fields to depend on color. This possibility arises in the instanton model of 
Carter and Diakonov jjj], who remarked that the Dyson-Gorkov equations close in color space on the condition that 
gapped and ungapped quarks can develop different masses. We studied the effects of this additional degree of freedom 
in the limit of zero chemical potential || and found no change in the phase structure for ratios B/A < N c /2. The 
situation is different for finite /i. Choosing different masses for the two gapped and the N c — 2 ungapped quarks leads 
in some limits to an increase in the pressure of phases with finite diquark fields. The main consequence is an increase 
of both the parameter range and the region in the (fj,, T) plane for which the mixed broken symmetry phase exists. 
This result is obvious since the mixed broken symmetry phase is the only phase which can exploit this additional 
degree of freedom. However, there is also a general mechanism at work by which a change of phase structure always 
takes place above a certain threshold for the coupling constants. 

We wish to illustrate these effects in a few cases and concentrate on the limit N c = 3 and m = for clarity. We 
denote the chiral fields for color 1 and 2 by a\ and that for color 3 by 03. In order to permit these fields to be 
different, we include the projection of the interaction onto a chiral- As channel as described in ||. For /i ^ 0, the 
thcrmodynamical potential becomes 



n(ai, <T 3 , A) - AA 2 + B{S5 1 a\ + ^g x g z + (3 3 a 2 3 ) - log[(ai - /1) 2 + A 2 + T 2 } - log^ + fi) 2 + A 2 + T 2 } 
-ilog[(a 3 -/i) 2 +T 2 ]-i 



1 log[(<7 3 - M) 2 + T 2 } - \ log[(<7 3 + fi) 2 + T 2 ], (41) 



where 



4 C 4 2C 1 C 

^9 + 3S' ^ = 9"3S' ft= 9 + 3S' (42) 

describe the coupling between a\ and cr 3 . C is a coupling constant associated with the chiral- As channel By fine 
tuning the various Lorentz and color channels which compose the random matrix interactions, we can realize a range 
of ratios B/C for any fixed B/A. New patterns of symmetry breaking only appears when C > and the chiral- As 
channel is attractive |jj . 



J This result is only true for a smooth cutoff F(q). For a sharp cutoff F(q) = <9(A — fi), where Q(x) is the Heavyside function, 
the diquark field exhibits the square root singularity. 
4 This range is -3/16 < B/C < 3/2 for B/A < 3/4 and then reduces linearly to B/C = 3/2 for B/A = 3/2. 
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Consider first single-gluon exchange. The coupling ratios are B/A = 3/4 and C/B = —3/16. The chiral- As channel 
is repulsive, and the phases with the largest pressure always satisfy <j\ — 03, which is the case shown in Fig. 2. If 
we now fine tune the interactions so as to increase B/C to B/C ~ 3/2 while keeping B/A = 3/4, the phase diagram 
changes substantially. The first-order line between the chiral and trivail phases splits at the point where it meets the 
diquark transition line into two first-order lines. Together with the T — axis, they delimitate a wedge of mixed 
broken symmetry phase with a\ ^ 173 and A ^ 0. Thus, with an attractive channel B/C ~ 3/2 the mixed broken 
symmetry phase appears much earlier than previously discussed. This is an extreme case. For fixed B/A, the mixed 
broken symmetry phase does not appear immediately as B/C increases from 0: there is a threshold value above which 
the new phase appears (this value is B/C ~ 1.29 for B/A = 0.75). Therefore, it takes large variations of the coupling 
constants B/A and B/C away from the values expected for single-gluon exchange to modify the phase structure of 
Fig. 2. 

To complete the picture of the effects of an attractive channel C > 0, consider next B/A = 3/2. This ratio can only 
be realized by a color diagonal interaction for which B/C = 3/2 ||. In this case, there is no freedom in fine tuning 
B/C to other values. The effect of splitting masses in colors is now maximal. (3% in Eq. ( [i"2"| ) vanishes and 03 decouples 
from the other two fields. Its gap equation leads to the same solutions as a pure chiral phase with a 7^ and A = 
and the phase diagram for the third color is that of Fig. 1. For B/A = B/C = 3/2, the partial pressure for colors 1 
and 2 has the same form as in the limit N c = 2, see Eq. (|37|). There is then a rotational symmetry between a\ and 
A for /1 = and o\ vanishes for \l 7^ while A = 2/ A — fi 2 — T 2 for /1 2 + T 2 < 2/ A. The overall phase diagram is 
obtained by superposing the two pictures; we have a mixed broken symmetry phase with <Ji — 0, 173 7^ 0, and A / 
on the left of the first-order line L Xt \ of Eq. (|l9|). This phase is contiguous to a pure diquark phase which develops 
as in Fig. 2 on the right of L x 1 and below the semicircle. With respect to Fig. 5, the mixed broken symmetry phase 
thus occupies a wider area in the (fx, T) plane. 

A color-6 condensate. It has been suggested that an energy gain may result if the third color condenses in a 
spin-1 color symmetric state ||. We find no such condensation for single-gluon exchange. We can however fine tune 
the interactions so as to keep B/A = 3/2 and allow this channel to develop. The Fierz constant for projecting on 
a color 6 is half that for a 3 color state. Thus, we can understand qualitatively how a color-6 behaves by repeating 
our previous analysis with an additional phase which now may develop inside a semicircle of radius \[2 smaller than 
the radius of the diquark phase. This semicircle crosses the first-order line beween chiral and diquark phases at 
B/A ~ 0.754 and a color-6 phase can, in favorable cases, increase the pressure of the system. However, as far as 
QCD is concerned, a color-6 condensate should be very small since its threshold ratio is very close to B/A = 3/4. 
This is an example of a result which cannot be regarded as robust: in a microscopic theory, the fate of the color-6 
phase will inevitably depend on the details of the interaction and on whether these give rise to exponential tails for 
the associated condensation field. By contrast, the chiral and diquark phases are fully developed for B/A = 3/4 and 
the phase diagram in Fig. 2 should thus be considered robust. 

V. CONCLUSIONS 

Our random matrix model leads to a thermodynamic potential which has a very simple form. Yet, it contains 
enough physics to illustrate in a clear way the interplay between chiral symmetry breaking and the formation of quark 
Cooper pairs. We have found that this interplay results in a variety of phase diagrams which can be characterized 
by a total of only six different topologies. Single-gluon exchange leads to the topology shown in Fig. 2, a phase 
diagram familiar from microscopic models. We have considered the chiral and scalar diquark channels, which seem 
the most promising ones. We have found that chiral and diquark phases are fully developped for the ratio realized 
by single-gluon exchange and that it takes large variations in the coupling ratios B/A and C/B to depart from that 
result. On the other hand, other less attractive condensation channels seem sensitive to coupling constant ratios and 
are expected in microscopic models to depend on the details of the interactions. Furthermore, these channels develop 
weak condensates at best. It is worth keeping in mind that the present picture is mean-field; we expect that quantum 
fluctuations will inevitably have large effects on the weak channels, which should thus not be considered robust. We 
conclude that QCD with two light flavors should realize the topology suggested by single-gluon exchange and that 
this topology is stable against variations in the detailed form of the microscopic interactions. 

Many effects lie in the mismatch between the number of colors involved in each order parameter: all N c colors 
contribute to the chiral field, while a Cooper pair involves two colors. A result of this for N c — 3 is that the chiral 
broken phase prevails at low densities. Furthermore, our model reproduces to a reasonable extent the expected 
limits at small and large N c . We expect such counting arguments to be valid in both microscopic models and lattice 
calculations. More generally, we believe that random matrix models can provide insight into calculations of QCD at 
finite density by providing simple illustrations for many of the mechanisms which are at work. 
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Not all the phase diagrams that we have studied are directly relevant to QCD. However, many of their characteristics 
such as the presence of tricritical and tetracritical points are generic to systems in which two forms of order compete. 
Our model could naturally be extended to the study of nuclear or condensed matter systems in which such competion 
takes place. The construction of a genuine theory seems technically involved at first glance M but in fact contains 
only three basic ingredients: the identification of the symmetries at play and their associated order parameters, the 
knowledge of the elementary excitations in a background of condensed fields, and the calculation of the range of 
coupling constants realized by the random matrix interactions. These three components are sufficient to build a 
thermodynamical potential and determine the resulting phase structures in parameter space. 
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APPENDIX A: CALCULATION OF ai(Nc) AND a 2 (N c ) 



In this appendix, we calculate the ratios ai(N c ) and 0:2 (-/V c ), which respectively characterize the onset of the 
^A-phase and the appearance of the tetracritical point. 



1. The ratio ai(JV c ). 



The mixed broken symmetry phase appears first for the coupling ratio B/A for which the instability line L m j x of 
Eq. (|25|) crosses the first-order line between chiral and diquark phases. This crossing occurs on the T — axis. The 
line L mix meets the T = axis at Eq. ©, 



2 

/^mix 



N„A N r 



AB 2 IB ' 



while the condition of equal pressure at that point gives 



il A -n x = (N c -2) + (A+ B)^ mbc - 2 log 



AUoj 



Bfi 



(Al) 



Combining these two equations gives the determining equation for a\ = B/A, 



N r ( N, 



4ai \2ai 



Nr lOE 



I f N r 



2 \2ai 



0. 



(A2) 



This is a transcendental relation which needs to be solved numerically for each N c . In particular, we find ot\{2) = 1, 
ai(3) = 1.05, and ai(N c 00) ~ 0.321iV c . 



2. The ratio a 2 (N c ). 



Two conditions determine a2(N c ). Coming from small B/A, a2(N c ) corresponds to the ratio at which the wedge 
of the xA-phase reaches the critical point (p^^Ti) of Eq. (p7|). Decreasing B/A from large values, a2(N c ) marks 
the coexistence of the tetracritical point (/i4,T4) and the tricritical point (At3m)T3 m ). We now show that these two 
conditions are equivalent, and thus that the transition from Fig. 4 to Fig. 5 is continuous. 

To proceed, we concentrate on the region near [n^T^) where we perform a small field expansion of the potential 
Q(a,A) inEq. (§), 

fl(a, A) m Q{0, 0) + a a 2 + b a A 2 + ^ a 4 + ^ A 4 + ci a 2 A 2 + G(min 6 (A, a)). (A3) 

Here, the coefficients of the quadratic terms are linear in 6 fi = fi — [14 and ST = T — T4 and vanish at (/X4, T4), while 
those of the quartic terms are 
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2B 2 

~n7 



A 2 N r . 



1? A " 



2AB A 2 

~Nr. 2~' 



(A4) 



to leading order in <5/i and ST. 

The condition that the xA-phase includes {n^^T^j can be determined as follows. While B/A increases from below 
&2(N C ), the tip of the wedge of the %A-phase slides along the the line of equal pressures between chiral and diquark 



phases, which we denote by L x _>a,i- Minimizing Q(a, A) in Eq. 
A-phases, this line is determined near (/Z4, T4) by 



= 0. 



to hnd the respective pressure in the x~ and 

(A5) 



The left boundary of the ^A-phase is the instability line 

d 2 n 



dA 2 



0. 



A=0,cr=er o 



where <7o is the a held in the chiral phase. From Eq. (A3), this gives near (ji^, I4) 

t ~ u a ° _ n 

iv m ix ~ Do 9C1 — U. 



(A6) 



(A7) 



The tip P m i x of the xA-phase is the intercept of L m i x and L x _»a,i- From Eqs. (A5) and (A7) P m i x thus obeys 

«o (c? - a 2 b 2 ) = 0. 



(A8) 



Now, imagine that B/A is strictly smaller but near ot2(N c ). Then, P m i X ^ (/X4, T4), and we have ao 7^ 0. The loca tion 
of Pmix is therefore determined by the condition a\b\ — cf, where the coefficients a\, b\ and c\ are those of Eq. (A4) 
augmented by linear corrections in Sfi and ST. If we now let B/A — > a2(N c ), the tip P m ix reaches (/Z4,T4) and the 
linear corrections in ax, &i, and ci vanish. Using Eq. (A4), and solving a\b\ = for P//1 gives 



, AT . B 4N C - n! /2 ^2N c - 4 

a2(iVc) = A = ^THq ' 

which is the result stated in Eq. (p8|). 

Coming from large ratios B/A, the determining condition for B/A = a2(N c ) requires that (^4^4) obeys Eq. (p^), 

^(a,A(a)) = 0, 



where A(tr) is the diquark held in the x^-phase. Taking the derivative of Eq. (A3) with respect to A, we find 
A 2 (er) pa — (ci<r 2 + b )/bf which when inserted back in Eq. (A3) gives 



d A VL 



4 



(a,A(a))wa? 



(A9) 



in the neighborhood of {ji^T^). Setting the right side to zero, we recover the previous condition of Eq. (AS). 
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FIG. 1. Phase diagram for B/j4 < 0.418 and jV c = 3, as a function of the quark chemical potential n and the temperature 
T. For this and all the following figures, continuous curves represent second-order lines while first-order lines are plotted with 
dots. Here, \ i s the chiral phase and A is the diquark phase. First- and second-order lines join tangentially at the tricritical 
point (n3,Ta). 
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FIG. 2. Phase diagram for the coupling ratio realized by single-gluon exchange, B/A — 0.75. The transition from the chiral 
phase x t° the diquark phase A is first-order while that from the A to the trivial phase is second-order. 
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FIG. 3. Phase diagram for a ratio B/A = 1.054. The mixed broken symmetry phase emerges out of the chiral phase 
via a second-order transition at /i m i x , and undergoes a first-order transition towards the diquark phase at jii. 
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FIG. 4. Phase diagram for a ratio B/A =1.1. Compared to Fig. 3, the tricritical point (fi3, T3) no longer exists. The trivial, 
chiral, and diquark phases now meet at the new critical point (^4^4). We note in particular that the first and second-order 
lines bordering the chiral phase do not meet tangentially at (/X4,T4). 
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FIG. 6. Phase diagram for a ratio B/A = 1.8. This topology is characteristic of large ratios B/A, which favor diquark over 
chiral condensation. The A-phase occupies a large part of the phase diagram, at the expense of the mixed broken symmetry 
phase. The chiral phase has completely vanished. We note that since Hermitian interactions realize ratios B/A < 1.5, this last 
topology actually describes a case out of reach for our model. 
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FIG. 7. Order parameters for the coupling ratio of single-gluon exchange, B/A — 0.75. The corresponding phase diagram is 
shown in Fig. 2. Chiral and diquark fields vanish continuously along the second-order lines. 
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FIG. 8. Expanded view of Fig. 3 near the wedge of mixed broken symmetry. We note that the second-order line between 
the x~ to the phase does not meet tangentially with the first-order line between the xA- and the A-phases, as indicated by 
the short dashed line. 
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FIG. 9. Expanded view of Fig. 3 near the tricritical point (/X3,T3), showing that there is a very short first-order segment 
between the chiral and trivial phases. 
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FIG. 11. Phase diagram for a ratio B/A = 0.2 and a small quark mass m. 
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FIG. 12. Phase diagram for the coupling ratio of single-gluon exchange, B/A — 0.75, and a small quark mass m ~ 10 MeV. 
The first-order line ends at a critical point (/n e ,T c ). 
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FIG. 13. Phase diagram for B/A = 1.0 and a small quark mass m ~ 10 MeV. 
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Table 1. A few characteristic coupling ratios. 



Ratio B/A 


N c = 2 


N c = 3 


iV c 


oo 


Onset of diquark condensation 




0.41s 


n 1 on 
0.139 


AT 

iv c 


C*l(iV c ) 


1.0 


1.050 


0.321 


N c 


Disappearance of the tricritical point. 




1.061 


0.354 


N c 


a 2 (N c ) 


1.0 


1.163 


0.354 


N c 


Maximum ratio N c /2 


1.0 


1.5 


0.5 


N c 


Single-gluon exchange 


1.0 


0.75 


0.5 
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